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Abstract 

Using 586 pb _1 of e + e~ collision data at Eqm = 4170 MeV, produced at the CESR collider and 
collected with the CLEO-c detector, we observe the process e + e~ — > ir + ir~ h c (lP). We measure 
its cross section to be 15.6 ± 2.3 ± 1.9 ± 3.0 pb, where the third error is due to the external 
uncertainty on the branching fraction of ip(2S) — > ir°h c (lP), which we use for normalization. We 
also find evidence for e + e~ — > rfh c (lP) at 4170 MeV at the 3a level, and see hints of a rise in the 
e+e" ->■ ir + TT~h c (lP) cross section at 4260 MeV. 
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In a previous Letter [1], the CLEO Collaboration investigated fifteen transitions to the 
J/if), ip(2S), and Xcj of charmonium states produced in e + e~ collisions with Ecu = 3970 — 
4260 MeV. The data were grouped into three energy bins (3970 - 4060, 4120 - 4200, and 
4260 MeV) roughly corresponding to the ?/>(4040), ■0(4160) and V(4260) regions. Increases 
in the e + e~ — > 7i + 7i^J/ip and e + e~ — > 7r°7r°J/'?/> cross sections at Eqm — 4260 MeV were 
attributed to V(4260) production j2j. In this Letter, we extend those investigations to 
search for 7r + 7r~, 7r°7r°, tt°, and 77 transitions to the h c (where h c = h c (lP)). We use the 
same 60 pb _1 of data with Eqm = 3970 — 4260 MeV (referred to as the "scan data" [3]) and 
the same energy binning, but we also now use 586 pb _1 of data collected at Eqm = 4170 MeV 
(referred to as the "4170 data"). The 4170 data set is an order of magnitude larger than 
was available at that energy for the previous study. The observation of transitions to the 
h c could provide insight into the perplexing nature of the charmonium states above DD 
threshold Q. It has also inspired new ways to search for and study bottomonium states, 
such as the hb (H 0] . 

We search for the processes e + e~ — > Xh c (X = 7r + 7r~, 7r°7r°, 7r°, 77) by reconstruct- 
ing the h c through 77/7 and the r\ c through: 2(7r + 7r~), 2(7r + 7r~)27r°, 3(n + n^), K ± K®ii Zf , 

K ± K° s 'K*'K+1t- , K+K-7T , K+K-7T+71-, K + K~7T + 7T~7T°, K + K-2{tT + 71-), 2{K + K~), 777T+7T-, 

and 7]2(tt + tt~), the same twelve modes used in the CLEO measurement of B(J/ip — > 7^ c ) 0- 
We also use a data sample of 24.5 million ip(2S) decays to reconstruct the process ip(2S) — > 
7T°h c using the same method. To eliminate dependence on the branching fractions of the 1] C , 
we take ratios of the cross sections (<7§) for e + e~ — > Xh c at center-of-mass energy E to the 
branching fraction (B§°) of ip(2S) -> ir°h c . We use B§° = (8.4 ± 1.3 ± 1.0) x 10~ 4 , measured 
by BESIII [8], to obtain erf. 

We utilize symmetric e + e~ collisions provided by the Cornell Electron Storage Ring (CESR) 
with center-of-mass energies at the ip{2S) mass and in the range 3970 — 4260 MeV. The 
resulting final state particles {K ± , 7r ± , and 7) are detected by the CLEO-c detector 0, 
which has a solid angle coverage of 93%. The momenta of charged particles are measured 
by concentric drift chambers [lOj, operating in a 1.0 T magnetic field along the beam axis, 
with relative momentum resolutions of ~0.6% at p — 1 GeV/c. To separate from tt^, 
two particle identification systems are used - one based on ionization energy loss (dE/dx) 
in the drift chamber and the other a ring- imaging Cherenkov (RICH) detector [ll[ . Photon 
energies are measured with a cesium iodide calorimeter, which has relative energy resolutions 
of 2.2% at £ 7 = 1 GeV and 5% at 100 MeV. 

We use standard track quality, particle identification, and calorimetry selection require- 
ments to reconstruct the exclusive processes e + e~ — > Xh c and ip{2S) — > n°h c with 
h c IVc- The 77's from the f] c are reconstructed in both their 77 and 7r + 7r~7r° decay modes, 
but the transition 77 from e + e" — > r]h c is only reconstructed in its 77 mode (due to large 
combinatoric backgrounds and small efficiencies for the 7r + 7r - 7r° mode). For ir° and 77 decays 
to 77, the mass of the pair of daughter photons is required to be within 3a of the nominal 
mass and is subsequently constrained to that mass. To reconstruct 77 — > 7r + 7r~7r°, the three 
pions must have an invariant mass within 30 MeV/c 2 of the nominal 77 mass. The Kg candi- 
dates are selected from pairs of oppositely charged and vertex-constrained tracks (assumed 
to be pions) with invariant mass within 15 MeV/c 2 of the Kg mass. In addition, we require 
that the photon from h c — > 777c cannot be paired with any other shower in the event to form 
a diphoton mass within 3a of the tt° mass. A four-constraint kinematic fit of all identified 
particles to the initial e + e~ four-momentum is then performed and the resulting fit qual- 
ity is required to satisfy xlc /d.o.f. < 5. This procedure sharpens the measured momenta 
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FIG. 1. The recoil mass of the 7r + 7r~ system versus the recoil mass of the jtt + tt~ system for 
candidate e + e~ — > ir + iT~h c ; h c — > ^r\ c events at Eqm = 4170 MeV. The signal appears at the 
intersection of the h c and r] c masses. The vertical lines indicate the region used to select the r\ c . 
The horizontal lines mark ±10 MeV/c 2 around the h c mass. 



in signal events and reduces backgrounds with missing or extra particles. For each decay 
mode of the 7] c , the candidate with the best fit quality is accepted. The selection criteria for 
ip(2S) — > 7i°h c are identical to that for e + e~ — > Xh c , except for an additional requirement 
suppressing ip{2S) — > 7i + 7i^J/ip by the exclusion of any event with a 7r + 7r~ pair with a recoil 
mass within 15 MeV/c 2 of M(jfy). 

We select the i] c by requiring the recoil mass of the 7X system be between 2930 and 
3030 MeV/c 2 . We then search for the h c in the recoil mass distribution of the X sys- 
tem. Figure [1] shows a histogram of the 7r + 7r~ and 77r + 7r~ recoil masses for the process 
e + e~ — > 7T + 7T~h c at Ecm = 4170 MeV. A clear accumulation of events can be seen near the 
intersection of the h c and 7] c masses, which marks the signal. Background from the initial 
state radiation process e + e~ — > r yijj{2S)] ip(2S) — > n + it~ J /ip appears as a vertical band at the 
J/ijj mass and is well-separated from the signal. Other backgrounds, studied with dedicated 
background Monte Carlo simulations, are smooth and are due to the light-quark continuum 
(e + e~ — y qq) or DD production, simulated with previously measured cross sections Q. 

The yield of e + e~ — > n + n^h c events at E C m = 4170 MeV is determined by fitting the 
7r + 7r~ recoil mass distribution, after selecting the 1] C , with two components. The signal shape 
is described by a double Gaussian with floating mass and normalization, but with widths 
fixed by signal Monte Carlo. The background shape is a freely floating first-order polynomial. 
The resulting fit is shown in Fig. E^a). We find 131 ± 15 signal events with a significance 
of more than 10a. The significance, here and in subsequent fits, is calculated from log- 
likelihood differences between fits with and without a signal component. The resulting mass 
from the fit is 3523.86 ± 0.48 MeV/c 2 (statistical errors only), which is 1.5 MeV/c 2 lower 
than the PDG 2010 value of 3525.42 ± 0.29 MeV/c 2 [12|. This discrepancy, however, is less 
than the uncertainty of the initial e + e~ energy 2 MeV) used in the kinematic fit, which 
directly affects the measured dipion recoil mass. 

Fits to e + e~ -> (7r 7r /tt° /r))h c at E CM = 4170 MeV and fits to e + e~ vr+rtc with 
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TABLE I. Yields (Ng), significances, relative fitting and shape systematic errors, efficiency ratios (R € ), normalized cross sections (cr^ /Bl), 
and cross sections (crj|) for each reaction e + e _ — > X/i c . The third error on cr^ is from B^f 0]. 
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FIG. 2. Fits to determine the yields of h c events from (a) e + e — > tt + tt h c at Eqm = 4170 MeV; 
(b) e+e" -> r/^ c at Eqm = 4170 MeV; (c) e+e" -> 7r + vr"/i c at E C m = 4260 MeV; and (d) the 
normalizing mode, ip(2S) — > 7r°/i c . 



-E"cm = 3970 — 4260 MeV follow the same procedure except that, due to lower statistics, 
the mass is fixed to the value obtained previously, 3523.86 MeV/c 2 . The resulting yields 
and significances are listed in Table [B We find >3<r evidence for e + e~ — > r]h c at 4170 MeV 
(Fig. |2](b)) and hints of a signal (2.6a) for e + e" 



-> 7i+7i-h c at 4260 MeV (Fig.^c)). 



The normalizing mode, ip{2S) — > n°h c , is also fit using the same method and with a 
floating mass (Fig. 121(d)). The yield is measured to be 202 ± 16 events. The resulting mass 
is 3525.27 ±0.17 MeV/c 2 (statistical errors only), consistent with, and highly correlated to, 
a previous measurement by CLEO using a similar method 13 . 



We calculate the ratios of the cross sections of e + e — > Xh c at energy E (<Xg) to the 
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branching fraction of ip(2S) — > TT°h c (£>^°) using: 

c e n£r: u 

where A^ is the number of ip{2S) decays, Ce is the luminosity at energy E, N~§ and N^f 
are measured yields, and R e is a ratio of selection efficiencies: that of e + e~ — > Xh c to that of 
ip(2S) — > 7i°h c . Since the ratio of efficiencies for each r) c decay mode is not perfectly constant 
(with 10% — 20% variations), we weight the individual efficiency ratios by the number of 
ip(2S) — > 7T°h c events we observe in each r] c decay mode, which we obtain through the fitting 
procedure described above. The errors on the efficiency ratios include errors due to Monte 
Carlo statistics and errors on these individual yields. 

Previously determined systematic errors are used for (2%) [14| and Ce (1%) [Hj]. 
Most systematic errors on individual track and photon reconstruction efficiencies cancel in 
the ratio of efficiencies, R e . However, for the transition particles, the X in the numerator 
and the 7r° in the denominator, a 1% relative error is assigned for each track and a 2% error 
for each photon. A conservative 5% systematic error is included for our determination of R e , 
which relies upon signal Monte Carlo distributed according to phase space. This systematic 
error is estimated by using extreme variations of the r\ c substructure - for example, by 
replacing 2(K+K~) by 0(1020)0(1020). 

Systematic errors in N§ and Nl° due to the fitting procedure are evaluated by varying the 
order of the background polynomials, varying the fit ranges, and varying the bin sizes. Based 
on Monte Carlo studies, we also use background shapes determined by xlc/d-o.f. sidebands 



'10 < xlc/d.o.f. < 35). For NI° , we alternatively use an ARGUS distribution 16] for the 
background. 

Systematic errors due to signal shapes are evaluated by varying the signal mass and width. 
The largest deviations occur when the signal widths are allowed to float. This variation 
determines the shape systematic error on Nl and N£ 17 q . For other A^, where the statistics 
are lower, the width variation is performed by scaling the width by the deviation observed 
between data and signal Monte Carlo in the fit for N£ 17 q , which is ~ 20%. Variations of 
the signal mass produce smaller deviations. 

The final numbers are listed in Table [B The TT + Tr~h c cross sections as a function of 
center-of-mass energy are summarized in Fig. [3J Notice that the 7i + 7i^h c cross sections are 
of a comparable size to those of 7t + tt~ J '/ip. There is also a suggestive rise in the cross section 
at 4260 MeV, which could be an indication of V(4260) production, but will require further 
data to be definitive. 

Projections of the n + TT~h c Dalitz plot at Eqm = 4170 MeV are shown in Fig. H]and are 
compared to phase space Monte Carlo. To separate signal from background, the number of 
signal iT + 7T~h c events in each bin is determined by the fitting procedure described above. 
The efficiency is relatively uniform across the Dalitz plot. More data would be required to 
investigate any possible discrepancies of the data with phase space. 

Assuming the E CM = 3970 - 4060 MeV and E CM = 4170 MeV data correspond to 
■0(4040) and -0(4160) production, respectively, we convert cross sections to upper limits on 
branching fractions using the same conversion factors listed in a previous CLEO analysis of 
this region [l|. The results are listed in Table [TTJ Assuming the 4260 MeV point is purely 
due to y(4260) production, we set a limit on its branching fraction to 7i + 7i^h c relative to 
7T+7T-J/0 of < 1.0 at 90% confidence level (C.L.). 
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FIG. 3. Cross sections as a function of center-of-mass energy. The triangle shows the cross section 
for e + e~ — > 7r + 7r~/i c at Eqm = 4170 MeV; the closed circles are for the same process at other 
center-of-mass energies. For reference, the e + e~ — > tt + tt~J/iP cross section [1] is indicated by open 
circles. The inner error bars are the statistical errors; the outer error bars are the quadratic sum 
of the statistical and systematic errors. 
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FIG. 4. The (a) tt + tt~ and (b) hcir^ mass distributions from e + e~ — > -ir + -ir~h c at Eqm = 4170 MeV. 
The points are obtained by fitting for the h c yields in bins of tt + tt~ or TT^h c mass. The histogram 
is signal MC, generated according to phase space and scaled by the total h c yield. 



In summary, we observe the process e + e~ — > -K + n~h c at Ecu = 4170 MeV and find its 
cross section to be comparable to the corresponding cross section for J/ip production. This 
has already resulted in new methods to search for and study the hf, using e + e~ collisions 
above BB threshold [6|. We also see hints of a rise in the 7r + TT~h c cross section at Eqm — 
4260 MeV. Further data will be required, however, to determine if this rise can be attributed 
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TABLE II. Upper limits (at 9 


0% C.L.) on branching 


fractions for the ^(4040) and V(4160) to Xh c . 
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to the ^(4260). 
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